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We systematically discuss properties of quantum disordered states of the quantum Hall bilayer 
at vt = 1. For one of them, so-called vortex metal state, we find ODLRO (off-diagonal long-range 
order) of algebraic kind, and derive its transport properties. It is shown that this state is relevant 
for the explanation of the "imperfect" superfluid behavior, and persistent intercorrelations, for large 
distances between layers, that were found in experiments. 
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I. INTRODUCTION 

Electrons in quantum Hall bilayer systems at total fill- 
ing factor ut = 1 naturally correlate in two different 
ways due to Pauli principle and Coulomb interaction. If 
the layers are sufficiently far apart, dominant correlations 
would be those of intralayer kind because electrons in one 
layer are unable to sense what is taking place in the op- 
posite layer. This does not hold, however, in the limit 
of small layer separation. Instead, with decreasing d/ls, 
the ratio of the distance between layers to the magnetic 
length, the correlations between electrons in different lay- 
ers gain strength and begin to compete with intralayer 
correlations. It is the interplay of those two kinds of cor- 
relations that we focus on in this paper. 

For the case of prevalent interlayer correlations, 
there are already a few theoretical models at hand 
which provide a satisfactory description: Ill-state given 
by Halperin's 1 111 wavefunction fm = n^o'C^T — 
z 3l) Uk<i( z kl~ z n) H m , n ( z mt ~ z nl), quantum Hall ferro- 
magnet?^ condensate of excitons^ or composite bosons. 4 
Nevertheless, both theoretically and experimentally, it is 
evident that with increasing d/ls a quantum disorder- 
ing of this state is bound to take place. For example, 
the tunneling peak observed by Spielman et al. 5 is in- 
deed sharp and pronounced, but its nature is more that 
of a resonance than of the speculated Josephson effect, 
while the temperature dependences of Hall and longi- 
tudinal resistance in experiments of Kellogg et al& and 
Tutuc et alZ do not provide support to the predicted 
Berezinskii-Kosterlitz-Thouless (BKT) scenario of a bi- 
layer finite temperature phase transition. 2 . Deeper under- 
standing of the regime d ~ Is is therefore an important, 
open problem in the physics of quantum Hall bilayers and 
strongly correlated electron systems in general. 

Hereinafter we present some results which pertain to 
quantum disordering that is believed to take place in 
the quantum Hall bilayer at vt = 1. The ground state 
at d — is a Bose condensate well-described by 111- 
wavefunction due to Halperin, while the low-lying excita- 
tions are composite bosons i.e. electrons dressed with one 
quantum of magnetic flux* 4 . The idea of disordering that 
we employ is to allow the formation of composite fermions 



(i.e. electrons dressed with two quanta of magnetic flux) 
that coexist with composite bosons^ There are two ways 
to introduce composite fermions into the Bose conden- 
sate and this will be explained in Sec.|TTJ We then pursue 
a phenomenological Chern-Simons transport theory a la 
Drude in order to examine the elementary predictions of 
those two model states. In Sec. IIIII we arrive at an ef- 
fective gauge theory for both cases. This enables us to 
calculate the correlation functions, modes of low-lying ex- 
citations and characteristic off-diagonal long-range order 
(ODLRO). We will be primarily interested in the pseu- 
dospin channel of these states. In one of those, the so- 
called vortex metal state that we believe may appear in 
the bilayer at larger d/ls as a manifestation of increas- 
ing intracorrelations, we derive an algebraic ODLRO. 
In Sec. IIVI we focus on the incompressible region and 
the crossover around the critical layer separation. We 
will argue that our field-theoretical, homogeneous pic- 
ture in fact suggests that vortex metal, if relevant for 
the strongly-coupled, incompressible region, may appear 
only localized, in form of islands in the background of 
the superfluid state for smaller d/ls- In Sec. IVl we give 
a more thorough analysis of the experiments on bilayer, 
addressing especially the compressible, weakly-coupled 
region, and the question of persistent intercorrelations- 
in the framework of the vortex metal state. Sec. I VII is de- 
voted to discussion and conclusion. For the sake of clarity 
and in order to make the text self-contained, some of the 
known results 8,10 will be re-derived in this paper. 

II. TRIAL WAVEFUNCTIONS FOR THE 
BILAYER 

Building on Laughlin's proposal for the wavefunction 
of a single quantum Hall layer,— the construction of 
Rezayi-Read wavefunction. 12 for v — 1/2 and Halperin's 
111-wavefunction for bilayer ^ we may formally imagine 
that there are two species of electrons in each layer (z, w) 
which are all mutually correlated through intracorrela- 
tions (within the same layer) and intercorrelations (be- 
tween opposite layers) , Fig. Q] 

Starting from the Ill-function of the Bose condensate, 
we will minimally deform it in order to include the com- 
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FIG. 1: Correlations between electrons in two layers. 

posite fermions. Given that each particle binds the same 
number of flux quanta and taking Pauli principle into 
account, this becomes a combinatorial problem with two 
solutions. In the first case: 

*i = ^in^T - zft) no*! - «u)n^t - ^ 
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The first line in this formula can be recognized as 
Ill-function, followed by two v = 1/2 separate layers 
($/'s denote the Slater determinants of free composite 
fermions), while the last two lines stem from the flux- 
particle constraint (all these correlations are depicted on 
the left hand side of Fig. [TJ. V,A denote projection to 
the lowest Landau level (LLL) and fcrmionic antisym- 
metrization (independently for each layer), respectively. 
In thermodynamic limit, the relation between the num- 
ber of particles and flux quanta reads^: 

N<& = N b1 + N bl + N n + Nfi 
= 2N n + N b] + N H = 2N fi + iV bT + N H . (2) 

N§, is the number of flux quanta through the system, and 
N blJ and Nf a are the number of bosons and fermions in- 
side the layer a respectively; a =t, I is the layer index. 
Eq. @ enforces an additional constraint Nff — Nf±. 
Therefore, the number of fermions is balanced in two lay- 
ers, while boson numbers are not subject to any such con- 
straint. This fact is important because of the broken sym- 
metry of spontaneous interlayer phase coherence in the 
Ill-state, which demands nonconservation of JV&j — iV&j. 
Although we work in a fixed (relative) number represen- 
tation (allowed in a broken symmetry case), to account 
for a broken symmetry situation we need to have a possi- 
bility of unconstrained relative number of bosons. Then a 



superposition of the wavefunctions of the form in Eq. ([T]) 
would lead to the usual representation. 

In the second case which is expected to describe dom- 
inant intracorrelations, fermions bind exclusively within 
the layer they belong to (right side of Fig. []} and the 
corresponding wavefunction is: 
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In this case, the flux-particle relationi^ is: 

Ng, = 2N f1 + 27V bT = 2N fl + 2N H 
= 2N f] + N b] -+N bl = 2N fl + N b] +N bl , (4) 

implying that both fermion and boson numbers must be 
balanced: N fl = N fl ,N b1 = N bl . 

In Ref. [1 the authors numerically calculated the over- 
lap of with the exact ground state wavefunction for 
a system of 5 electrons in each layer with varying d/ls- 
Their results seem to demonstrate convincingly that (at 
least for small systems) the approach with trial wave- 
functions that interpolate between two well-established 
limits, namely those of Ill-state and decoupled v = 1/2 
layers, is not only an artificial mathematical construction 
but corresponds to physical reality. Despite the fact that 
the number of electrons in this simulation is certainly 
well below the thermodynamic limit, the fact that the 
overlaps between and the exact ground state display 
peaks very close to 1 at small d/ls provides confidence in 
the choice of wavefunction (at least for small d/ls)- 

If there is a phase separation in between the sea of com- 
posite bosons and composite fermions, the phase transi- 
tion will be of the first order. Such a scenario is launched 
in Ref. 1 13L where the authors imagine static, isolated re- 
gions of incoherent phase inside 111 phase. Although 
this model correctly explains some of the observed phe- 
nomena (e.g. semicircle law), the persistence of inter- 
correlations in the weakly-coupled, compressible regime 9 
which gradually die out suggests a continuous transition. 
Such a possibility is naturally present in the picture of 
composite boson-composite fermion mixture. 

A transport theory of Drude kind can be easily formu- 
lated 8 if we consider that composite fermions bind two 
quanta of magnetic flux, unlike composite bosons which 
bind only one quantum of magnetic flux. As long as we 
remain in the RPA approximation, they can all be treated 
as free particles moving in the presence of the effective 
field which is given by the sum of the external and self- 
consistently induced electric field. In the first case (4'i), 
the effective field as seen by particles in the layer a is: 

£] =E°~2eJZ-e(Jl+Jt), (5) 
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where J^(m denote Fermi- and Bose-currents in the layer 
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Transport equations are: 
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f(b) - Pf{b) J f(b) 



(7) 



and, as required by symmetry, Pf^) = P 2 f{by whhe the 
total current is given by 3 a = 3% + Jy . We define single 
layer resistance (p 11 ) and drag resistance (p D ) as follows:: 

E^^J 1 , (8) 

(9) 



E 2 = p D j\ 



When both layers have the same filling, v\ = v% = 1/2, 
tensors pb,Pf are diagonal (because the composite par- 
ticles are in zero net field): p b = diag [p bxx , pbxx], Pf = 
diag [pfxx-i Pfxx] , and in the case of drag we have in ad- 
dition J 2 = 0, J^finitc. Then from Eqs. ©-© via ele- 
mentary algebraic manipulations we obtain: 
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Formulas in Eqs. (fT2"|) - (fT4")) include parameters S, p bxx 
and Pfxx, the last two being free parameters about which 
nothing can be said a priori. This prompted Simon et 
al£ to reason as follows. At large d/ls the number of 
composite bosons is small because the condensate is bro- 
ken and pbxx is large compared to 8, which is the typical 
Hall resistance. On the other hand, from the experi- 
ments^ we know that for large d/ls holds pf xx <C 8. 
Furthermore, even as oI/Ib is decreased, we expect pf xx 
to increase only slightly 4 All in all, for large oI/Ib they 
assume pbxx 3> <5 3> Pfxx, and if in addition we allow 
PbxxPfxx <C 5 2 , asymptotically we obtain: 
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Semicircle law follows directly from the previous formu- 
las: 



(Pxx) + (Pxy ~ o) 



4/, 



(18) 



in agreement with Ref. [lj (semicircle law is of general 
validity for twocomponent systems in two dimensions and 
it serves us as a crucial test for the the line of reasoning 
quoted above, which may at first sound somewhat naive). 

In the opposite limit (when d/ls is reduced), p bxx <C 
Pfxx "C S, because p bxx drops as a result of Bose conden- 
sation^. When pbxx —> we obtain the quantization of 
Coulomb drag: 
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as measured by Kellogg et aL— 

Let us return now to the case of dominant intracorre- 
lations, the vortex metal stated represented by Eq. ([3]). 
From Fig. Q] the formulas for effective fields arc modified 
into: 



£J = E CT -2eJ^-2eJ b r , 
=E° -e(Jl+J 2 + 2r f ), 
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and the analogous calculation yields the resistivity ten- 
sors: 
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In this case as well, there are two physically significant 
limits depending on the assumptions for the values of 
Pbxx, Pfxx- In the case when p hxx < Pf xx < S: 
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and the semicircle law follows, Eq. (fl8|) , whereas \p xx \ = 
|. Similarly, in the regime pi, xx <C S <C p/ xx we de- 



duce the quantization of Coulomb drag, Eqs. (fT9| - ([20ll 

We emphasize that these two limits arc different from 
those in Simon et a/4 For example, the semicircle law 
was derived assuming that pb xx is small (which is ex- 
actly the opposite situation to the one in Ref. [1), while 
Pf xx is not necessarily small with respect to 5. As noted 
in the first case above, the exact values for Pb xx ,Pf xx 
are in fact unknown and this prevents us from discrim- 
inating between the different proposed limits. In other 
words, we cannot say which one of the proposed limits 
is plausible - the analysis above serves us only to con- 
clude that each of the two composite boson-composite 
fermion mixed states are able (with certain assumptions) 
to reproduce the phenomenology of drag experiments. 



III. CHERN-SIMONS THEORY FOR BILAYER 

Encouraged by the preliminary analysis from the pre- 
vious section, we will pursue the idea of composite boson- 
composite fermion mixture further by formulating an ex- 
ample of Chern-Simons (CS) field theory which can con- 
tain wavefunctions ^1,^2 a s ground states. We do not 
embark on such a task only for the sake of completeness, 
but also because such a theory would enable efficient cal- 
culation of response functions and provide insight into 
the long range order of the system and the nature of low- 
lying excitations. General drawback of CS theories is the 
inability to include the projection to LLL which is the 
arena where all the physics must be taking place. Never- 
theless, we will use these theories established in the works 
of Zhang et al^ of composite bosons and Halperin et 
al*£ for composite fermions because even analyses done 
in advanced, projected to the LLL type of theories, in 
the work of Murthy and Shankar— came to the conclu- 
sion that in order to get, in the most efficient way, to the 
qualitative picture of the physics of response, the usual 
CS theories are quite enough and accurate. In addition 
to this simplification, in constructing the CS theories, we 
will neglect the antisymmetrization requirement implied 
by Eqs. |T]) and ([3]). The reason for this is that just like in 
hierarchical constructions, c omp osite fermions represent 
meron excitations, see Ref. Ildl . that quantum disorder 
the Ill-state and, as it is usual when we discuss the dual 
picture of the FQHE 18 , we do not extend the antisym- 
metrization requirement to the quasiparticle part of the 
electron fluid. 

Therefore we start from the lagrangian given byi^: 
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where a enumerates the layers, \5 a i $ CT are compos- 
ite fermion and composite boson fields in the layer a, 

v n = V U = V «.i V U = V 11 = V e, and t he densi- 
ties are 8p a = Sp F + 5p B . By A (and B) here we 
mean external fields in addition to the vector potential of 
the uniform magnetic field, Ab, which is accounted for 
and included in gauge fields a^K Therefore we have 
a F(B)a = k F(B)<y _ Ab Externa i ne id s A and B cou- 
ple with charge and pseudospin, and in general we must 
introduce four gauge fields a F ( B ) cr . Fortunately, not all 
of them are independent. In the first case, the relation 
analogous to Eq. ((2|) becomes the following gauge field 
equation: 
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From the equations above, it is obvious that there 
are only two linearly independent gauge fields: &c — 
a + a = - — J*— and a s = - — ^— , and Eq. (32|) ex- 
pressed in Coulomb gauge reads: ^^f- = 8p^ + Spi = Sp 
and = 5p F ^ — 8p F ^ = Sp F (ac, as are the trans- 
verse components of the gauge fields). These are the 
constraints we wish to include into the functional inte- 
gral via Lagrange multipliers a^ and Og . The interaction 
part of the lagrangian is easily diagonalized by introduc- 
ingVc = Y^ andVs = ^. 

The strategy for integrating out the bosonic functions 
is the Madelung ansatz (p a = p a + ~p a e l9 < J , which ex- 
pands the wavefunction in terms of a product of its am- 
plitude and phase factor, while fcrmionic functions are 
treated as elaborated in Ref. [l6|. After Fourier trans- 
formation, within the quadratic (RPA) approximation, 
and introducing substitutions Sp c = Sp^ + Sp^,Sp s = 



5 pi — Sp^ ,i = F,B and 
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^, all the 
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terms neatly decouple into a charge and a pseudospin 
channel: 
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where 8a$ = a$ — A , Sac = ac ~ A, Sa^ = Og — B , 
Sag = as — B, p b -mean density of bosons in (each) 
layer. In writing down Eqs. ([33 ]) . (|34|) we utilize a com- 
pact notation suppressing k (—k) dependence, where 
all the quadratic terms of the type [X + Y] stand for 
[X(-fc) + Y(-k)] [X{k) + Y(k)]. K 00 (k) and K n (k) are 
the free fermion (RPA) density-density and current- 
current correlation functions r& In the long wavelength 
limit (k/kf -C 1) they can be explicitly evaluated from 
the general expressions^: 
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where x = jrp;, fc/-Fermi wavevector and ny-fermion 
density, 0- Heaviside step function. The mass appear- 
ing in expressions for -Koch Kn is equal to the bare elec- 
tron mass only in RPA approximation (in which we work 
here) . 

Focusing on the charge channel only, Eq. (f3"3")l . and 
integrating out first Spc, then Og and Sac, we arrive at 
the density-density correlator: 
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then k 



i—t-x, and we conclude that as Li 



ix),Ku w 
-► (and 

0) the system is incompressible in the charge 
channel, so long as there is a thermodynamically signifi- 
cant density of bosons ~p b . 

In the pseudospin channel, we are primarily looking 
for the signature of a Bose condensate i.e. whether there 
exists a Goldstone mode of broken symmetry and what is 
the long range order of the state. Therefore, in Eq. 
we set = = 0, and integrate over a$ , as, Sps ■ 



(0s(-k)9 s (k)) 



V s 



»2 jl 
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where a — jKqq — Kn( 2 ^-) 2 (in Appendix [Al we give the 
full linear response in the pseudospin channel). Indeed, 
there exists a Goldstone mode, albeit with a small dissi- 
pative term (which, if desired, can be removed by pairing 
construction^): 
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Even for large x, it is easy to check that the pole re- 
mains at the same value if we assume ~p b 3> nt (which 
is, in fact, the most appropriate assumption in this 
case). Also, the imaginary term disappears in this case. 
Such robust Goldstone mode implies the existence of a 
true ODLRO and the genuine Bose condensate. Gold- 
stone mode Lu°(k), Eq. (|39|) . is easily observed in Fig. [2j 
where we plotted the real part of density-density cor- 
relation function 7Too(fc), Eq. (|A1[) . in terms of parame- 
ters Q = k/kf and x = u)/{kkt). Other (fixed) param- 
eters are: m = Ib = 1, d = 0.5, e = 12.6, V$ = ird/e, 
Pb + n f = 1/(4tt), n = n f /p b = 1/10. 
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FIG. 2: Re7roo(fc) and the Goldstone mode in the case of 

Let us return to the second case, that of Eq. Q and 
dominant intracorrelations. According to Fig. [TJ rela- 
tions Eq. ([32]) are modified to become: 
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It is obvious that in this case we have only 3 linearly inde- 
pendent gauge fields, namely ac = - — — = - — ^ — , 

- . Introducing the same 



as 



and aps = — 



2 ■= — 2 

substitutions as before, the lagrangian again decouples 
into a charge channel: 
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where 5oq S 



Oq — Bo,SaFS = &fs — B, all the other 
symbols have retained their meanings. 

This time we will not analyze the charge channel in de- 
tail. To this end, we note that the system in incompress- 
ible in this sector, the fact which is easily established by 
integrating out all the gauge fields, densities and boson 
phase in Eq. (HT]) . 

In the pseudospin channel, a calculation of the density- 
density correlator leads to the conclusion that in this 
channel the system is compressible (see also Fig. [3]). The 
9 — 9 correlator is: 
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2K n . For 



small k/kf and x, the correlator diverges for uj° — -^Eh. = 
const, which obviously contradicts the original assump- 
tion for the range of x and hence we reject this pole. For 
x >■ 1 (and still k <C fc/), the relations Eq. (|3"5|) . are 
approximately K aa « —^s,Ku ~ and we obtain 

two poles: 



(44) 



n\k) = ^^J\+r ] + \jlT±n, (45) 

where n — ~p b /nf is the ratio of boson to fermion density 
(Eqs. (|4"4"|) . (|4"5"1) hold for any n, although in the physical 
limit that we are presently interested in, n may be re- 
garded as small). In Fig. [3] we plotted the real part of 
the density-density correlation function in the case of 5*2, 
Eq. (|A4[) . In contrast to Fig. [51 here we opt for ui and Q 
as free parameters and set d = 1.5 and n = "p^/nj = 1/10 
as the more likely values in this case. Distinctive feature 
of Fig. [3] at ijj = 1 is the plasma frequency fl° and the 
smaller singularity at u = 1/10 corresponds to w°. There 
is also a striking absence of Goldstone mode in this case. 

We now proceed to calculate ODLRO in the pseu- 
dospin channel of ^2- As it turns out, ODLRO will 
be nontrivially modified and assume algebraic form. We 
know that interaction does not affect the value of char- 
acteristic exponent^ and therefore set Vs = 0. Bearing 
in mind that we work in the long wavelength limit, we 
arrive at the following expression for the correlator: 
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where we introduced top — 7Trif . After contour integra- 
tion over u*i^: 
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ODLRO: 



„ \ , , which leads to the algebraic 



r'|/(n) 




FIG. 3: Re7roo(fc) in the case of ^ 



„/|-(l-2r / +o(r / 2 )) 



(47) 



This algebraic ODLRO persists as long as 77 > (func- 
tion / is positive everywhere in this domain). The ex- 
pression Eq. (|47j) is formally reminiscent of BKT XY- 
ordering, only the role of the temperature is overtaken 
by the parameter n (the analysis of this paper assumes 
temperature T = 0). Pursuing this analogy further, we 
conclude that the relative fluctuations of composite bo- 
son and composite fermion density represent the mecha- 
nism which may lead to the ultimate breakdown of the 
111-condensatc. 



IV. EVOLUTION OF THE GROUND STATE 
WITH d 

In order to investigate the transition from the incom- 
pressible, Ill-like state at lower d/l B , to the compress- 
ible, possibly vortex metal- like state at higher d/l B , we 
are motivated to introduce what we call generalized vor- 
tex metal. In addition to the ordinary vortex metal (^2), 
we include (in each layer) another kind of composite 
fermions that connect to the composite boson sea as in 
the case of The generalized vortex metal is clearly 
the only additional option left of connecting electrons di- 
vided in composite bosons and composite fermions beside 
the two extreme cases, and 4 , 2- Again, in this state 
some of composite fermions connect in the manner of in- 
state to the composite bosons and the rest of composite 
fermions connect exclusively to the composite bosons of 
the same layer in the manner of the Rezayi-Read state. 
This is succinctly represented by the following gauge field 
constraints: 



y x a Ba - A - ■■ ■'- ^<>''-' l> -'- '>a7 2) 



1 

2^ 



5 PB1 + S PBL + SpW + Sp™ + 26p [ ;l,(48) 
Vxflf^ 5 PB] + 6 P Bi + 2Sp { pl + 2^^,(49) 
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—v x 4' 

Z7T 



25p B a + 28p%l + 2^(50) 



where the superscripts (1), (2) indicate composite 
fcrmion species in each layer. Chern-Simons theory easily 
follows from the above gauge field equations and yields 
incompressible behavior in the charge channel. In the 
pseudospin channel: 



C P8 = K&\5a$Zr + K£(8agM) 



FS\2 



-(2), 



,FS\2 



m m 
Zn Zn 



(2)/ 



,FS\2 



(51) 



where the linearly independent fields are given by as — 

n Bf_ n Bl „ F T_„-F4. „„ „FT _ 



— , subscripts 
1,2 distinguish between composite fermion species and 
s denotes antisymmetric combination of the densities in 
two layers (like in Sec. IIII[) . A noteworthy feature of the 
lagrangian, Eq. (|5Tj) . is the existence of Vhc, the hard- 
core repulsion term between the two species of compos- 
ite fermions inside each layer. The presence of such a 
term (added by hand) is natural if we imagine composite 
fermions residing in two separate Fermi spheres. How- 
ever, the danger of blindly introducing this term is that 
it may incidentally bring about the incompressible be- 
havior (otherwise not present) in the system. We have 
verified that this is not the case here i.e. the system 
remains incompressible whether or not we choose to in- 
troduce Vhc- It therefore appears more intuitive to keep 
Vhc, taking the limit Vhc — > oo in the end. Step by step, 
eliminating all the gauge fields, we are lead to the follow- 
ing correlation function: 



(0s(-k)9 s (k)) = 



(52) 



and the low-energy spectrum is dominated by the plasma 
frequency: 



An 



u (k) = — y / p b n f2 , 



(53) 



where n f2 is the density of the composite fermions which 
bind exclusively within the layer they belong to. Gen- 
eralized vortex metal therefore is a state that only sup- 
ports gapped collective excitations, despite the presence 
of composite bosons and the kind of composite fermions 
which enforce interlayer correlation. If it is pertinent 
to the region of the tunneling experiments of Spielman 
et al£ and counterflow experiments of Kellogg et al&, 
we believe that our homogeneous theory of Sec. IIIIIIVI 



then suggests that (generalized) vortex metal can appear 
only as localized islands (due to presence of disorder at 
low temperatures) amidst the background of phase 
(Fig. rj|. In Fig. 2] depicted are weakly-coupled vortex- 
antivortex pairs i.e. meron-antimeron pairs (due to the 
charge degree of freedom, there are four kinds of merons 2 ) 
inside the vortex metal phase. They are expected to exist 
in the vortex-metal phase on the grounds of disordering 
of the correlated phase. As argued in Ref. |To|, the inclu- 
sion of composite fermions into the 111 state (\I/i and vf^) 
corresponds to creation of meron-antimeron pairs. There 
are more pairs and more of larger size as d increases con- 
sistent with the BKT picture of the phase that supports 
algebraic ODLRO, Eq. gZ}. 




d < d c 



d> d c 



FIG. 4: Evolution of the ground state with varying d, before 
and after the transition at d — dc- The regions with meron 
pairs represent the vortex metal (^2) phase. The background 
represents the superfluid (^1) phase. 



V. FURTHER COMPARISON WITH 
EXPERIMENTS 

In this section we wish to address in depth the potential 
of the model states, Vl/i and in explaining the phe- 
nomenology of experiments on bilayer. The key question 
in this analysis is: what is the nature of the compress- 
ible phase corresponding to higher d/ls that still harbors 
some of the intercorrelation present at lower d/ls 9 ? 

The answer to this question cannot be given by look- 
ing at simple transport properties. In Sec. |TT] it was 
shown that both and ^2 in certain regimes can re- 
cover the two main experimental findings of Kellogg et al. 
in drag experiments: the semicircle law 9 and the quan- 
tization of Hall drag resistance^. On the other hand, 
our Chern-Simons RPA approach at T = stresses that 
all that states considered in this paper are incompress- 
ible. However, at finite T, a finite energy^ is needed 
to excite a meron in ^ and therefore ^2 seems like a 
better candidate for exhibiting compressible behavior at 
any finite T or, at least, a very small gap. Furthermore, 
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within the vortex metal picture, ^2 allows the follow- 
ing simple scenario. For p bxx <C Pf xx "C S, one gets the 
semicircle law derived in Sec. HU As d/ls increases, the 
density of bosons decreases and one enters the regime 
Pfxx -C Pbxx <C S where \p xx \ > p xy (as witnessed in 
the experiments^). The persistence of enhanced longi- 
tudinal drag resistance^ up to very high d/lg provides 
additional support to our choice of ^2 which can explain 
the remaining intercorrelation (drag) in the case where 
explicit tunneling is absent. Finally, as pb xx — > 00, both 
resistances go to zero, the bilayer decouples and bosons 
vanish from the system. 

Our picture is certainly incomplete because it does not 
explicitly include the effects of disorder (which must be 
very relevant for the physics of bilayer in the regimes 
d ~ Ib - a simple way to see this is to look at the be- 
havior of measured counterflow resistances^ p xx , p xy 
that enter the insulating regime very quickly after pass- 
ing through vt = 1)- Fertig and Murthy2i provided a 
realistic model for the effects of disorder and in their 
disorder-induced coherence network in the incompress- 
ible phase of the bilayer, merons are able to sweep by 
hopping across the system, causing the activated behav- 
ior of resistance (dissipation) in counterflow. This finding 
is consistent with our own. 

At the end our picture is in the spirit of the Stern and 
Halperin proposal^ but instead of the 1 /2 compressible 
phase coexisting in a phase separated picture with the 
superfluid phase we assume the existence of the 

vortex metal phase (^2)- This coincides with the Fertig 
and Murthy proposal^! for the incompressible region that 
explains the "imperfect" superfluid behavior, ft is the 
continuous extrapolation of this phase separated picture 
that brings and favors ^2 for larger d/ls (instead of 
There "I^ is able to explain the persistence of intercorre- 
lations through enhanced longitudinal drag accompanied 
by the absence of tunneling and phase coherence^ 

Finally, we are able to account for the effects of the 
layer density imbalance in tunneling, drag^ and counter- 
flow^ experiments. Spiclman et al. 22 observed that small 
density imbalance stabilizes the resonant tunneling peak 
- a simple reason for this is that \&i can easily accom- 
modate the fluctuations in density (see comment after 
Eq.©). Because of the same reason, Hall drag resistance 
remains quantized up to larger d/ls in the presence of 
density imbalance. On the other hand, the enhancement 
of longitudinal drag resistance at large d/ls was also re- 
ported^ to be insensitive to density imbalance. While the 
reason for this cannot be seen only from looking at the 
form of ^2 (this state constrains both fermion and boson 
numbers in two layers, see comment after Eq. [J) , we be- 
lieve that meron excitations are responsible for absorbing 
the density fluctuations, especially at finite T. 

Recently, the quantum Hall bilayer was probed using 
resonant Rayleigh scattering^ for samples with differ- 
ent tunneling amplitudes and when the in-plane mag- 
netic field is present. They detected a nonuniform spa- 
tial structure in the vicinity of the transition, suggesting 



a phase-separated version of the ground state. Our re- 
sults (for zero tunneling limit and excluding disorder) 
hint that such phase separation may indeed be neces- 
sary to invoke in order to achieve a full description of 
the strongly-coupled, incompressible phase and the tran- 
sition in a bilayer. 



VI. DISCUSSION AND CONCLUSION 

In conclusion we showed how two model states, ^1 
and \&2j can account for the basic phenomenology of the 
bilayer that came up from various experiments. 

A very interesting question pertains the model state 
^2- As effectively a state that represents a collection of 
meron excitations interacting through topological inter- 
actions, a question comes when they are in a confined 
(dipole) phase and when in a metallic (plasma) phase. 
So in principle we can expect that the static correla- 
tor in Eq. (|4"7)) can be reproduced by considering a 2D 
bosonic model with meron excitations interacting via 2D 
Coulomb plasma interaction^ Therefore we believe that 
the Laughlin ansata^ of considering (static) ground state 
correlators as statistical models in 2D can be applied 
here also. We expect that the ground state correlators 
in a dual approach, in which we switch from compos- 
ite fermion to moron coordinates, can be mapped to a 
partition function of a 2D Coulomb plasma^ The 2D 
Coulomb plasma has two different phases. For large (3 
(inverse T) the charges form dipoles and the system is 
with long range correlations (no mass gap). At some 
critical j3, dissociation of dipoles occurs and we have a 
plasma phase with a Debye screening, and therefore a 
mass gap. Thus calculations that will capture more of 
the meron contribution than our RPA approach in the 
Chern-Simons theory may find a transition and expo- 
nential decay of the correlator (Eq. (|47p). before reaching 
~p b = limit. Indeed our ODLRO exponent in Eq. (|47)) 
at ~p b = is 1 which is well above the exponent of 
the BKT transition or critical exponent 1/4. At that 
point our system may develop a gap in the pseudospin 
channel and completely lose interlayer coherence (expo- 
nential decay of correlators). Furthermore, we expect 
that the superfluid portion of the composite boson den- 
sity will disappear leading to compressible behavior in 
the charge channel^ This is all consistent with experi- 
ments^^ which find that, at d/ls ~ 1.84, the vanishing 
of the conventional quantum Hall effect and the system's 
Josephson-like tunneling characteristics occurs simulta- 
neously. Intercorrelated bosons continue to exist without 
a superfluid property and lead to enhanced v = 1 drag 
at large d/ls- They disappear from the system around 
d/l B w 2.62. 
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APPENDIX A 

In order to extract the response functions in functional 
integral formalism, one needs to integrate over all degrees 
of freedom except those of the external fields. The inte- 
gration of these fields in the RPA approximation prover- 
bially reduces to the Gaussian integral: 



d(z, z*)exp [—2 * wz + u * z + vz* 



-exp 



, Rew > 0. 



For the pseudospin channel in the case of (Eq. (TT|)) 
we therefore obtain the following linear response: 



1 



7T o(fc) - ^(0(3^ s + 2a) + l), 



(Al) 



1 — Z7T 

n i{k)=n w (k) = ^—K 11 {l + nVs), (A2) 

7m(fc) = l{2V s (i - nv s )(K u - Bh.) 

a m 

+K u /K 0Q -4a&}, (A3) 
m 



where a = - K n {^) 2 , fi 



V s 



and 



4-"-00 "iiv, k 

E = V s (l - QVs) + 2a. 

The response functions in the case of the pseudospin 
channel of * 2 (Eq. are: 

-oo = (A4) 
1 ik 

7T01 = 7T10 = X^~ A ' ( A5 ) 



tt, , = i{A 2 + A(2K n - ^ - WW\^r)}„ i A(, ) 



, 2np b 



2A' n n v 2t ^ 



A EE 



where VK 4 ee 

- 2^ n + V s (±r, A ee W 4 (^) 2 - 2K n . 
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